Angular Dependence of the Superconducting Transition Temperature in 
Ferromagnet-Superconductor-Ferromagnet Trilayers 
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The superconducting transition temperature, T c , of a ferromagnet (F) - superconductor (S) - fer- 
romagnet trilayer depends on the mutual orientation of the magnetic moments of the F layers. This 
effect has been previously observed in F/S/F systems as a T c difference between parallel and antipar- 
allel configurations of the F layers. Here we report measurements of T c in CuNi/Nb/CuNi trilayers 
as a function of the angle between the magnetic moments of the CuNi ferromagnets. The observed 
angular dependence of T c is in qualitative agreement with a F/S proximity theory that accounts for 
the odd triplet component of the condensate predicted to arise for non-collinear orientation of the 
magnetic moments of the F layers. 

PACS numbers: 74.45.+c,75.70.-i,74.78.Fk 



Electronic properties of a superconducting (S) film can 
be significantly altered by direct contact with a metallic 
ferromagnet (F) 0, 0] . Penetration of the superconduct- 
ing condensate into a ferromagnet results in the breaking 
of singlet Cooper pairs by the exchange field and leads to 
a decrease of the superconducting transition temperature 
(T c ). Recent discoveries of quantum interference effects 
in F/S heterostructures have shown that the F/S prox- 
imity effect is more interesting and complex than just the 
reduction of T c 0-0. Examples of non-trivial F/S prox- 
imity phenomena include (i) oscillations of T c with the F 
layer thickness (df) in F/S multilayers 0, 0|, (ii) oscilla- 
tions of the critical current in S/F/S Josephson junctions 
with d f @ and (in) manipulation of T c in F/S/F @, Q 
and F/N/F/S [Ii] structures (N denotes a nonmagnetic 
metal) by switching the magnetic state of the F layers. 

The superconducting condensate of a F/S system is 
predicted to consist of a singlet component and a triplet 
component with unusual symmetry (Tol - [l6j |. The triplet 
part of the condensate function is expected to be s-wave 
(i.e. even in the momentum variable). Therefore, by 
the Pauli exclusion principle, it must be odd in time. 
The s-wave character of the condensate makes it robust 
against impurity scattering, and thus the odd triplet pair- 
ing should survive even in dirty F/S systems. If the mag- 
netization of the ferromagnet is not spatially uniform, the 
triplet part of the condensate has three non-zero compo- 
nents corresponding to three projections, S z = (0, ±1) of 
the Cooper pair spin. The S z = ±1 components are im- 
mune to pair breaking by the exchange field, and thus can 
penetrate deep into the ferromagnet 17]. Recent mea- 



surements of the critical current in Josephson junctions 
with noncollinear magnetic barriers support the existence 
of this long-range odd triplet superconductivity 
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The simplest model system for studies of the effect of 
spatially non-uniform magnetization on superconductiv- 
ity is a F/S/F trilayer with non-collinear (neither paral- 



lel nor antiparallel) orientation of the magnetic moments 
of the F layers 0, EH- However, experimental tests of 
the F/S proximity effect in this system have been lim- 
ited to the collinear geometry [til Il9l - [22| . In this Letter, 
we describe experimental and theoretical studies of T c in 
CuNi/ Nb/ CuNi trilayers as a function of the angle be- 
tween the CuNi layer magnetizations. We compare our 
experimental results to theories of the F/S proximity ef- 
fect that take into account the odd triplet component of 
the condensate, and find that the clean limit of the prox- 
imity effect theory gives good qualitative description of 
the data. 

The multilayer samples for our experiments were grown 
by dc magnetron sputtering in a vacuum system with a 
base pressure of5.0xl0 -9 Torr. We deposited a series of 
multilayers of the following composition: thermally oxi- 
dized Si substrate/ Ni 80 Fe 2 o(4 nm)/ CuNi(d/)/ Nb(d s = 
18 nm)/ CuNi(d/)/ Ni 80 Fe 20 (4 nm)/ Ir2 5 Mn 75 (10 nm)/ 
Al(4 nm) with df ranging from 2 nm to 18 nm (see 
Fig.QJa)). The ferromagnetic Ni 80 Fe2o layers placed next 
to the CuNi layers ensure uniform magnetization of the 
CuNi layers in our measurements Q. The magnetiza- 
tion direction of the top F layer is pinned in the plane 
of the sample by exchange bias from the antifcrromag- 
netic Ir25Mn75 layer. The CuNi alloy was grown by co- 
sputtering Cu and Ni targets, and the Ni concentration 
of 52 atomic percent was calculated from the Cu and Ni 
deposition rates. 

We pattern the samples into 200 /im-wide Hall bars 
and perform four-point magnetoresistance measurements 
in a continuous flow 4 He cryostat. Fig.QJb) shows resis- 
tance versus temperature for the parallel (P) and an- 
tiparallel (AP) configurations of the CuNi layers for the 
sample with 5 nm thick CuNi. The superconducting tran- 
sition temperature difference between the P and AP con- 
figurations, AT C , is 7 mK. Fig. QJc) shows the resistance 
as a function of magnetic field applied along the exchange 
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FIG. 1. (color online), (a) Schematic of the sample. Here 
2a is the in-plane angle between the magnetic moments of 
the CuNi layers, (b)-(d) Experimental data for the sample 
with the CuNi layer thickness of 5 nm. (b) Resistance versus 
temperature for parallel (P) and antiparallel (AP) alignments 
of the CuNi layers, (c) Resistance versus in-plane magnetic 
field applied parallel to the pinned layer magnetization at T= 
2.745 K. (d) Resistance versus 100 Oe in-plane field angle, 2a, 
measured at the same temperature. 



bias direction in the middle of the superconducting tran- 
sition at T= 2.745 K. This figure demonstrates that a 
magnetic field of ±100 Oe is sufficient to reverse and 
saturate the magnetization of the free layer and thereby 
create well-defined P and AP configurations. Fig. [ljd) 
shows resistance versus the angle between the F layer 
magnetizations, 2a, measured at T= 2.745 K. For this 
measurement, we stabilize the temperature in the mid- 
dle of the superconducting transition to within ±0.1mK, 
rotate the 100 Oe external magnetic field in the plane 
of the sample and measure the sample resistance as a 
function of 2a. 

The angular dependence of T c is calculated from the 
angular dependence of the resistance, R(2a), such as that 
shown in Fig. [TJd) 23] . In order to convert the resistance 
data into T c , we use the slope of the R{T) curve, (4^), 
in the middle of the superconducting transition (at T= 
2.745 K for the sample with 5 nm thick CuNi shown 
in Fig. Hfb)). The angular variation of T c , AT c (d f , 2a) 
— T c (df,2a) — T c (df,0), is calculated for a given df as 
AT c (d fl 2a) = (R(2a)-R(0))/(j§). Fig. Ufa) shows 
AT c (df,2a) for samples with CuNi thickness dt from 2 
nm to 10 nm. This plot demonstrates that AT c (df,2a) 
significantly deviates from a simple dependence given by 
sin 2 (a). This is quantified in Fig. EJb), which shows the 
percent deviation of the normalized angular variation of 
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FIG. 2. (color online), (a) Angular variation of T c defined 
as AT c (df,2a) = T c (d f ,2a) - T c (d/,0) for samples with the 
CuNi layer thickness df from 2 to 10 nm. Inset: normalized 
AT c (df,2a) for df — 2 nm (squares) and sin 2 (a) (line), (b) 
Percent deviation of the normalized angular variation of T c , 
AT c (d f ,2a)/AT c {d f , 180°), from sin 2 (a). 



We next discuss the comparison of our results with 
theory. The transition temperature T c (df,2a) can be 
calculated for our system, which is in the clean limit, us- 
ing methods that include lH, 16| the odd triplet correla- 
tions that are generated by the presence of magnets with 
non-collinear magnetizations. These methods are based 
on self-consistent solution of the microscopic Bogoliubov- 
de Gennes 2J| equations. The matrix elements of these 
equations with respect to a suitable basis containing a 
sufficiently large number N of elements, can be evalu- 
ated analytically as a function of d s , a, and the other 
parameters in the problem. The eigenvalue problem for 
the resulting matrix is then numerically solved, using 



[16l ] a self consistent iterative process. To find T c , the 
self-consistency equation can be linearized [25| near the 
transition, leading to the form Aj = ^2 q Ji q A q , where the 
Aj are expansion coefficients of the position dependent 
pair potential 25] in the chosen basis and the Ji q are the 
appropriate matrix elements with respect to the same ba- 
sis, as obtained from the linearization procedure. These 
matrix elements can be written as Ji q = (Jf q + Jf q )/2, 
where, 
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cant for samples with df < 4 nm. The magnitude of the 
deviation reaches 20% for the df = 2 nm sample 
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FIG. 3. (color online). Experimental data (crosses) and theo- 
retical fittings (circles) of T c in the P state as a function of the 
CuNi layer thickness df. The inset is the comparison of AT C . 
Note the different scales for experimental data (left scale) and 
theoretical fittings (right scale). 



Here 7 = 7o/(47rD), with 70 being the usual di- 
mensionless superconducting coupling, e^"-'' are un- 
perturbed particle (hole) energies, D is the total 
sample thickness d in units of the Fermi wave- 
length, and No denotes that the sum is cut off 
at energies beyond the Debye frequency. We have 



,,,,,,, n v 2^^^ r K.qp r (dl n pC^ lr + ^mp c nr): fcqpr 

(2/cf) 3 / 2 Jq dzf(z) sin(fc g z) sm(k p z) sin(fc r z), where k q = 
qir/d, the z axis is normal to the layers (see Fig. [TJ and 
f(z) is unity in the superconductor and zero elsewhere. 
The quantity e± denotes the kinetic energy due to mo- 
tion in the transverse direction. The cfl and c?°l are the 



expansion coefficients of the quasiparticle amplitudes 25] 
in terms of the basis set. All energies are scaled by Ep. 
All the sums and integrals in the above equations can be 
evaluated numerically with great precision as a function 
of temperature and of the other parameters in the prob- 
lem. The value of T c is then found (25H27I as the highest 
temperature for which the largest eigenvalue of the ma- 
trix Ji q is unity. This largest eigenvalue is readily found 
numerically. 

In evaluating T c (df, 2a) we have used the actual value 
of d s (18 nm), and T c o = 4.5 K, where T c q is the bulk 
transition temperature of the material forming the super- 
conducting layer. The only other parameters in the prob- 
lem are [16|,|25| the superconducting coherence length £ , 
a parameter I characterizing the magnet strength (1=1 
in the half metallic limit and 1 = for a non-magnetic 
metal) and the dimensionless barrier strength (28|, Hb, 
characterizing interfacial scattering between the S and F 
layers. The intricacies of the numerics and the computing 
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FIG. 4. (color online). Comparison between the experiment 
(crosses, left scale) and the theory (circles, right scale) of the 
angular dependence of AT c (d/, 2a) for samples with different 
CuNi layer thicknesses df. For all df, experimental and the- 
oretical results have similar shapes despite the mismatch in 
absolute values. 



time required do not allow £o, and Hb to be used as 
free parameters in a best fit procedure. We have instead 
explored for best values within an experimentally reason- 
able region of parameter values. Results presented here 
are for £o = 16 nm, reasonable for a Nb film, / = 0.032 
appropriate for a relatively weak magnet, and Hb = 0.75, 
a barrier of moderate strength. 

The results are shown in the next two figures. In the 
main plot of Fig. [3] we show the results for T c (df,0), 
the transition temperature at a = 0, versus df. The 
experimental results show the expected [25| oscillations 
with df, which are reproduced quite accurately and in 
detail by the theory. In the inset, we show the differ- 
ence AT c (d/, 180°) between the transition temperature 
at 2a = 180° and at 2a = 0°. This quantity is also 
plotted as a function of df. Here, the general shape of 
the theoretical curve is qualitatively similar to the ex- 
periment, but the experimental overall scale of the phe- 
nomenon is much smaller than the theoretical one. We 
will return to this point below. 

We consider next the dependence on 2a. In the six 
panels of Fig. [4] we plot experimental and theoretical re- 
sults for AT c (df, 2a) as a function of 2a for six values of 
df in the experimentally relevant range. The shapes of 
the curves are similar to the experiment, although some 
deviations are present for thin (df < 4 nm) CuNi films. 
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However, the discrepancy in the magnitude found above 
at 2a = 180° affects, unavoidably, all the other angles. 

Thus, the agreement between theory and experiment 
is in all aspects satisfactory except for the overall scale of 
the T c shifts, which is much too large in the theory. This 
is a puzzle, since for apparently much more difficult prop- 
erties such as the shape of the curves in Fig. |4] and the 
main plot in Fig. U2 theory and experiment are in good 
agreement. One would think, a priori, that while the 
good detailed agreement in these quantities does not pre- 
clude discrepancies elsewhere, such discrepancies could 
not possibly be so large. We have also attempted to fit 
the data to the dirty limit case, using the methods of 
Ref. IT3I . In the dirty limit, the solution is mathemati- 
cally much simpler and we are able to perform a best 
fit to the data using six fitting parameters. Despite ob- 
taining in some cases unreasonable parameter values, the 
agreement with experiment for T c (df) (the experimental 
data shown in the main plot of Fig. [3]) is worse and so 
is the shape of the curves in the inset of Fig. [3] and in 
Fig. 01 particularly at smaller values of df. Even then, 
the discrepancy as to the size of the effect remains just 
as large and in some cases even larger than in the clean 
limit calculation. So, the discrepancy can not be due to 
the samples being closer to the dirty limit than thought. 
This large discrepancy in the magnitude of AT C between 
theoretical predictions and experiment has been noted 
before [1, 0, 0] in other contexts. 

It is natural to hypothesize that much larger surface 
scattering values, Hb, and smaller coherence length, £ , 
leading to a diminished proximity effect, would drasti- 
cally reduce AT C . This is indeed so: we find that increas- 
ing H B and decreasing £ sufficiently leads to values of 
AT C of the right order of magnitude, but then the depen- 
dence of the results on df and a is incorrect. A possible 
explanation is that something in the sample interfaces 
weakens the normal and Andreev scattering processes re- 
sponsible for the proximity effect in such a way that only 
the overall scale of the amplitude of the phenomenon is 
changed, but we cannot speculate on what this could be. 

In conclusion, we made measurements of the super- 
conducting transition temperature T c in CuNi/Nb/CuNi 
trilayers as a function of the angle between the CuNi 
magnetic moments and the thickness of the CuNi layers. 
We found that T c is a monotonically increasing function 
of the angle between the magnetic moments of the CuNi 
layers. For thin (df < 4 nm) CuNi layers, this function 
significantly deviates from a simple trigonometric func- 
tion sin 2 (a). The dependence of T c on the CuNi thick- 
ness and the functional form of the angular dependence 
of T c are reasonably well described by the proximity ef- 
fect theory in the clean limit but the overall scale of the 
angular variation of T c given by the theory exceed the 
experimentally observed values by two orders of magni- 
tude. Our measurements of the angular dependence of 
T c provide new experimental constraints on theories of 
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